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Functional Renormalization Group approach to the Kraichnan model 
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We study the anomalous scaling of the structure functions of a scalar field advected by a random 
Gaussian velocity field, the Kraichnan model, by means of Functional Renormalization Group tech¬ 
niques. We analyze the symmetries of the model and derive the leading correction to the structure 
functions considering the renormalization of composite operators and applying the operator product 
expansion. 


I. INTRODUCTION 

Understanding the scaling of the velocity structure 
functions in turbulent flows is a challenging problem. 
Such scaling is known to be universal in the inertial range 
and shows small departures from the behaviour predicted 
by Kolmogorov [T] . To make progress in understand¬ 
ing the anomalous scaling of the structure functions a 
model for an advected passive scalar field, the Kraichnan 
model [2], turned out to be illuminating (similar models 
were considered long time ago by Corssin and Obukhov 
mmy One of the features which makes this model ap¬ 
pealing is that the corrections to the canonical scaling 
of the structure functions can be analytically computed 
considering the zero modes of some operator Ad* which 
acts on the equal time correlation functions of the scalar 
field [U [B]. This result has been obtained via various 
methods among which the expansion in the Holder expo¬ 
nent appearing in the covariance of the velocity held (i.e.: 
the £ appearing in equation an expansion for large 
spatial dimension and the Batchelor limit (e close to 2) 
[M3- Furthermore the zero modes have been understood 
as statistical conservation laws [TO]. For a review of the 
results and a complete discussion we refer to m- In this 
work we study the scaling of the structure functions of 
the Kraichnan model by means of the Functional Renor¬ 
malization Group which offers the advantage of being a 
very adaptable framework. 

The Renormalization Groups proved to be ideal to un¬ 
derstand the anomalous scaling in statistical systems but 
the applications of its inherent concepts go far beyond. 
The Functional Renormalization Group (FRG) is a func¬ 
tional realization of the Wilsonian renormalization pro¬ 
gram and is applied successfully to many problems re¬ 
garding out of equilibrium systems, for a review see [12] 
and m- In the FRG framework one considers a scale 
dependent effective action, the Effective Average Action 
(EAA), which interpolates between the classical action 
and the standard effective action mm- Within this 
framework the issue of a fixed point for the functional 
integral associated to the randomly stirred Navier-Stokes 
equation has been first considered in m and more re¬ 
cently in [TE] and [HI [2D]. In particular in |2U] a fixed 


point solution of the FRG equation associated to this 
problem has been investigated in d = 2 and d = 3. Nev¬ 
ertheless there is not still a quantitative prediction for 
the anomalous exponent of the structure functions and 
it may turn useful to address this question in a simpler 
setting. Let us stress that the FRG has been proved to 
be a solid framework in similar problems such as the scal¬ 
ing exponents of the Kardar-Parisi-Zhang equation and 
its correlation functions in the strong coupling regime 
[211 122] (see [23] for the Burgers’ equation). 

The RG analysis of the Kraichnan model has been 
performed in a perturbative framework by means of an 
e —expansion which is very similar to the one used in the 
analysis of critical phenomena |24l . The £—expansion 
framework has been applied also to many other problems 
ranging from the randomly stirred Navier-Stokes equa¬ 
tion to magnetic hydrodynamics, we refer to [231 and [26] 
for a complete list of references. Actually this framework 
turned out to be very efficient for the Kraichnan model 
and the anomalous exponents of the structure function 
have been computed up to order O (e 3 ) m\- The inverse 
and the Wilsonian RG and their connection to the zero 
modes results have been studied in [28; . In this work we 
analyze the FRG approach to the Kraichnan model; this 
allows us to investigate a number of issues in a simpler 
setting and shows that the FRG is a viable framework 
for the computation of the anomalous exponents for the 
structure functions. Thus the main purpose of this work 
is to study in a concrete example the techniques required 
to understand the scaling of the structure functions in 
the FRG framework and in particular the role of com¬ 
posite operators. The paper is organized as follows: in 
the section [TT] we introduce the path integral formalism 
for the model. In section m we consider the relevant 
RG equations and show how to obtain the latter ones 
within the FRG approach. In section m we study the 
IR fixed point of the model while in section [V] we con¬ 
sider the anomalous dimension of the relevant composite 
operators. In section [Vi] we put together our results and 
discuss the scaling of the structure functions. In the final 
section we summarize our findings and consider possible 
outlooks. 
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II. PATH INTEGRAL FORMULATION OF THE 
KRAICHNAN MODEL 

In this section we recall the definition of the Kraichnan 
model and derive its path integral representation. The 
symmetries of the model are described in appendix [A] 
together with some non-renormalization theorems which 
follow. Let us consider a scalar quantity 0 (t, x), e.g.: the 
temperature, which is advected by a fluid having random 
velocity v (i, x). The equation for the dynamics of 9 (t, x) 
is 

d t e + v i d i 9-^d 2 e = f (i) 

where n is the molecular diffusivity and / is a stochas¬ 
tic forcing term with Gaussian statistic and zero average. 
Note that equation ([l]) is a stochastic differential equa¬ 
tion with multiplicative noise and is interpreted in the 
Stratonovich sense. More precisely the covariance of the 
forcing term is: 


(/ / (i 2 , x 2 )) = S (ti - t 2 ) F 

— Dg 


We have: 

Z = j VfVvV9P ff P vv 5 (d t d + v'cM) 
det (dt + v l d l - |<9 2 ) 



where the determinant comes from the change of vari¬ 
ables in the path integral. For the time being we shall 
neglect this determinant and we will come back to it in 
a moment: 


Z = J VfVvV9e~^ fUf) 1 f e ~^S v{ ' Vv) - 


5 (d t 9 + v%0 - ^d 2 9 - / 


= j V9VfVvV9e~^ i f{ff) Vg-HW) 
i J ~9(d t 9 + Pd l 9-^d 2 9- /) 
= f V9VvV9e + ^ iS{n)S e-^ iv( ' vv) ~ lv 


exp 


'(d t < 


- I 9{d t 9 + Fd l 9- -d 2 t 


where Lp is the characteristic scale of the forcing (even¬ 
tually we will denote LJ . 1 = M). The random velocity 
v (i, x) is also Gaussian and has zero average: 


Vi(t 1 ,x 1 )v j (t 2 ,x 2 )) = d S(t 1 t 2 ) . 

(27r) J 

[ d d ke ikxi - ikx2 

P 

r ij 

(2) 

(fc 2 + m 2 ) 5+ff 

= D v 

(3) 


where 


Pi, 


— 5ij 


kikj 

~W 


is the projector onto the transversal components (the 
fluid is incompressible). Throughout this work we as¬ 
sume the Stratonovich discretization prescription for all 
the relavant formulas unless otherwise stated. Note that 
the mass m in the denominator avoids IR divergences and 
is a further scale besides Lp, which appears in the force 
correlator. Actually in the inertial regime the structure 
functions do not depend on m and we set m to zero (if 
one wanted to it is easy to keep this parameter). We will 
not need to introduce m at all since the FRG framework 
is safe from IR divergences by construction. 

Averages can be computed via the following path in¬ 
tegral: 


In the last line we redefined the field 9 via 9 —> id. The 
determinant in equation Q is very similar to the one 
encountered in the path integral representation of the 
Langevin equation, see for instance [25] and references 
therein. Such determinant depends on the convention 
chosen for Heaviside theta function in zero: if 9 (0) = 0 
we are considering the Ito calculus while if 0(0) = 1/2 
we are considering Stratonovich. In particular, as shown 
in | i28| . using Ito discretization prescription the determi¬ 
nant does not enter in the path integral and the action 
has the same form of the action Sg defined in equation 
([5]) below but with a different molecular diffusivity. As 
we already said in this work we will follow Stratonovich 
discretization prescription and we shall include the deter¬ 
minant contribution in the action via two Grassmann odd 
fields which we refer to as ghosts. Such choice may seem 
somehow more complicated since it doubles the num¬ 
ber of fields. Nevertheless we would like to keep such 
ghost fields since they carry a physical and geometrical 
meaning in the path integral approach to classical me¬ 
chanics mm and make manifest the supersymmetry 
in stochastic processes, see [32] and references therein. 
The ghost action associated to the determinant is sim¬ 
ply: 

S 9 h = J c (d t + v l dt - c . 

The full action is thus composed of two parts: S 9 h and 


Z = 


j VfVvV9P f f P vv S (0- 9 0 ) 



\o D §0 + \vD v h 


15) 


where 9 0 is a solution of equation (TT|) and Pf f and P vv are 
the Gaussian weighting factors of / and v respectively. 


Note that the field 0 is the response field present in the 
Martin-Siggia-Rose formalism [33], see also [MIES], and 
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that Sg and S g h are very similar due to the fact that 
equation ([I]) is linear in 9. Thus the stochastic problem 
in equation 0 has been mapped to the field theory com¬ 
ing from expression Now we are able to exploit the 
techniques available for standard field theory and in par¬ 
ticular we will employ the Functional Renormalization 
Group discussed in the next section. Finally let us recall 
the dimension of the various quantities entering in the ac¬ 
tion, we have [9] = T 1 / 2 , [0] = L~ d T~ 1 / 2 , [A,] = L~ e 

and [k] = kp. Following [24] we define g = D 0 /n which 
has dimensions [< 7 ] = L~ e . 


III. FUNCTIONAL RENORMALIZATION 
GROUP FOR THE KRAICHNAN MODEL 


Functional Renormalization Group provides a non per¬ 
turbative framework for the implementation of the Wilso¬ 
nian integration of high momentum modes. For a general 
introduction we refer the reader to the following list of 
reviews [37lI41| . We will work via the Effective Average 
Action (EAA) which is a scale dependent generalization 
of the standard effective action. One first define a modi¬ 
fied generating functional of connected Green’s function 
W k : 

e W k [J] = j Vxe -S-AS k +J X 

where A S k suppresses the integration of momentum 
modes p 2 < k 2 and is quadratic in the fields with a kernel 
Rk, i.e.: A S k = \ J xRkX- Note that this acts like an 
infrared cutoff. Thus let T k be the Legendre transform 
of W k and define the EAA subtracting the cutoff action 
which we added at the beginning: 

Tfe = Ffe — A Sk ■ 


The k— dependence of the functional T k satisfy the fol¬ 
lowing exact equation [141416] 




( 6 ) 


where dt = kd k is the logarithmic derivative with respect 
to the cutoff and is the Hessian of the EAA. To 
concretely employ equation © we will need to resort to 
some approximations and we will make an ansatz for T k . 
Such procedure has been proved robust in many fields, 
especially to determine the scaling of statistical system 
at criticality, see ^ 2 ] for an overview. 

Equation | 6 |) will provide the running for the couplings 
in our ansatz. Nevertheless we are also interested in the 
running of composite operators. To derive the flow equa¬ 
tion for composite operators let us consider the one loop 
case. Besides the cutoff term let us add to the action a 
source e coupled to a composite operator O , then the one 
loop EAA reads 


r M = S + eO + ^Tr log + R k + eO< 2 >) 


where S^ and O ^ are the Hessians for the action and 
the composite operator respectively. If we now derive the 
above expression with respect to the scale k we obtain: 

d t r M = ±Tr (s< 2 > +R k + SOM)' 1 d t R k . 

Finally, in order to derive the running of the composite 
operator, we just need to take a functional derivative with 
respect to the source e and set this to zero: 

d t O = ~Tr (V 2 ) + A) -1 OM (S^+R.y 1 6 t R k . 

The exact flow equation for the running of the compos¬ 
ite operator can be found from the above expression just 
by “RG-improving” the one-loop case (exactly as it hap¬ 
pens with the EAA). To see this explicitly we consider 
that, given a modified action S + eO, both the connected 
Green’s function generating functional W k [J, e] and its 
Legendre transform [p,e\ depend on the source e. In 
particular from the definition of T k and T k we have 

r fc [<p, e] = J ipJ - W k [J, e] - A S k [y>] . 

In order to extract information regarding the composite 
operator O we simply take a functional derivative with 
respect to e and observe that 


S 

5e (x) 


r k W,e\ 


£=0 


SW k [J,e] 
Se (x) 


£—0 


This equation tells us that the running composite oper¬ 
ator [0} k is directly related to the Tfc [ip, e] via a func¬ 
tional derivative with respect to e. The flow equation for 
Tfc [p,e\ is simply given by equation ([ 6 ]) where the usual 
EAA is replaced with the e —dependent EAA. Thus the 
running of the first functional derivative of [ip, e] with 
respect to e can be directly obtained from ([ 6 ]) and, setting 
e = 0 at the end, one obtains: 


d t O k 






•(7) 


Equation 0 will give back a differential equation to 
which a suitable boundary condition is associated. Typ¬ 
ically one defines the composite operator of interest at 
some scale A. Clearly the RG evolution generally pro¬ 
duces mixing among different operators and the compos¬ 
ite operator has to be parametrized choosing a suitable 
ansatz. For a detailed discussion on the flow equation for 
composite operators we refer to [4Tj . 

Let us come back to the Kraichnan model. We choose 
the cutoff kernel R k as 


Rk 


( 0 R kM 0 0 0 \ 

R k gg 0 0 0 0 

0 0 R k , vv 0 0 

0 0 0 0 R kiC c 

V 0 0 0 R krcc 0 / 


( 8 ) 
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The matrix acts on the following vector in field space: 
(<?, 0, v, c, c). As common in the application of FRG to 
non-equilibrium systems we have chosen an off-diagonal 
cutoff for the response field and its companion (e.g.: 0 
and 9 respectively). Moreover we do not cutoff the fre¬ 
quency since we will be able to perform integration over 
frequency analytically. Note that adding the cutoff may 
break some of the symmetries analysed in appendix [A] 
Of course at the fixed point the scaling behaviour does 
not depend on the form of the cutoff but breaking some 
symmetries in the middle of the flow means that more 
terms have to be taken into account. It is clear that if 
we consider 

Rk = ^-(Af , (9) 

e^> -1 

where A = — 8 2 , none of the symmetry is broken and in 
particular shift symmetry is preserved. If we chose the 
optimized cutoff [43] 

R k = ( k 2a - A Q ) 9 (k 2 - A) (10) 


the renormalization procedure so that only finite terms 
remain. At this point we note that only operators like 
8989 are relevant in (11), i.e.: operators which are not 
multiplied by overall factors of T. In full analogy one 
can argue that for S 2n the relevant composite operator is 
(8989) n . This reasoning is at best heuristic and can be 
better justified by diagrammatic considerations as pre¬ 
sented in [24] where it is shown that only composite op¬ 
erators with at most 2 n 9 —fields appear in the OPE of 
S' 2 n ■ Here our main message is that the OPE suggests 
that the scaling of the structure functions is encoded in 
the scaling of some particular composite operators whose 
leading term in a derivative expansion is (8989) 2n . Note 
that this is the operator you get from a simple Taylor 
expansion of Stra¬ 
in order to derive the formula for the anomalous di¬ 
mension of a composite operator O we consider the 
Callan-Symanzik equation and implement renormaliza¬ 
tion multipling the operator by a coupling Zo■ Let us 
consider the coupling appearing in Sg and S g , the Callan- 
Symanzik equation reads: 


shift symmetry would be broken. 

We are interested in the scaling of structure functions 
S n of the field 9 at equal time, i.e.: 

S 2n = ((9(t,x)-9(t,y)) 2n ). 

Clearly S 2n is given by a product of several different com¬ 
posite operators and to determine the scaling behaviour 
of S 2n is a rather difficult task. It is convenient to rea¬ 
son in terms of single composite operators expressing Sbn 
via the Operator Product Expansion (OPE), or Short 
Distance Expansion. Such program works nicely in the 
e —expansion framework in j 24] and we will follow closely 
this procedure. Here the main idea is to identify the op¬ 
erators which determine the scaling in such a way to be 
able to apply the FRG framework afterwards. First of 
all let us consider the simplest case of S 2 and recall that 
the canonical dimension of the field is [6] = T 1 / 2 . For 
dimensional reason thus we have: 


(ndfj, + /3 k 8 k + 70 ) (O) = 0 (12) 

where we used the fact Dq is not renormalized and thus 
does not enter in the equation and 7 o = d t Zo IZq ■ More¬ 
over from dimensional analysis (in space and time) we 
have 


(■ ~D 0 8 Do - k8 k + T8 t + d£) (O) = 0 (13) 

+ (2 - e) D 0 8 d 0 

+ 2 k8 k + xd x — M8 m + do) P) = 9 (14) 


where d,Q and cIq are the frequency and mass dimensions 
of the operator O. Now let us plug (12) in (14), we obtain 


(find* + 7o + (2 — e) Dq8d 0 

+2 k8 k + xd x — M8m T do') P*) = 9 • 


Now we insert (13) in the above relation and find 


S 2 = T ■ Ci (htT, fix) + C 2 (htT, fix) [9 2 ] 

+T~ 1 C 3 (htT, hx) [0 4 ] -I- 1 - Di ( htT , fix) x 2 [8989] 


(/3 k ,8 k + 70 - eD 0 8 Do + 2 T8 t 

+ 2 dp + xd x — M8m H~ do) (O) = 0 . (15) 


+T 1 D 2 (htT, fix) x A 


(8989Y 


( 11 ) 


Note that odd terms in 9 are not allowed by 
Z 2 —symmetry. Moreover both the path integral and the 
correlation function in S n are shift invariant. As a conse¬ 
quence the operators entering in the OPE are also shift 
invariant, i.e.: only derivatives of 9 appear and all the 
terms in the first line, except the first one, vanish. Fi¬ 
nally let us recall that we are considering equal time cor¬ 
relators and so T must be set to zero. Now, if we make 
the hypothesis that the time dependence of the various 
terms entering in is not dramatically different from 
the canonical one, we can observe that for T —> 0 some 
terms vanish, some others stay finite and some diverge. 
We assume that the divergent terms are eliminated from 


Finally it is convenient to re-express the flow via the cou¬ 
pling g = Dq/k and to observe that we can rewrite the 
beta function of k via g exploiting the fact that Dq is not 
renormalized: 


/3 K d K = — k8 k 

K 


(*) 


hJL [ K ' 
^ dk 1 " 


D 0 

= - idgHd K . 

9 


k-^r (g x ) 

k8 k = dkK \ ’ k8 k 


Now we replace k8 k using equation (131 and we express 
f3 g via its dimensionless counterpart g = gk~ e . Note that 
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(3 g = (eg + (3g) k e and so 


—/3 s a id K — ~z (eg + (3g ) (—Dq8d 0 + T8t + dq>) 
= e (—DqDdo + TOt + do) 


at a non-Gaussian fixed point. Now we plug this in equa¬ 
tion (fl5|): 


(—e (—D 0 dr> 0 + T8 t + do) + 7 o ~ £-Do<5d 0 + 2T9t 

+2do + xd x — M8m + do) (O) = 0, 

(( 2-e)Td T + (2-e)d% 

-\-xd x — M8 m + do + 70 ) (O) = 0 . (16) 


This equation will be used in section [VI] to derive the 
scaling of various composite operators and to determine 
the scaling of the structure functions. Summarizing we 
have introduced the FRG framework for the couplings 
and for the composite operators. Thus we used the OPE 
to identify which is the most important operator which 
contributes to the scaling of the structure function. 


IV. THE IR FIXED POINT 

In this section we derive the IR fixed point whose scal¬ 
ing will be studied in the next sections. The presence 
of an IR fixed point for the Kraichnan model is known 
from e— expansion in [21] . In [21| it is argued that the 
beta function of g , which is the coupling we choose to 
parametrize the RG flow, is one loop exact. Thus we will 
first perform a one loop computation within the FRG 
framework and afterwards we will comment on how the 
fact that the beta function of g is one loop exact can be 
understood from an FRG perspective. 

We start by considering the ansatz: 



(191 


has the form 


G k = 


( Ggg Ggg 

Ggg 0 
0 
0 

V 0 


0 
0 

0 G vv 


0 

0 

0 

0 

~Gr 


0 \ 
0 

0 


G c 


0 / 


In particular we have 

Ggg = G, 


ee 


1 


-~ 2 ^ R k 


Gw — -; 

(<? 2 r 

_ 

C.C. 'GX nn 


iu + f q z + 
Do 


_ p 


R k 


ito + ^q 2 + R k 


Schematically the flow equation for two point function 
has the following form: 


r fc = 


in 


i 


i 


d t + v l di — —d 2 ) 9 — —QDgO + -vD L v 


+c(dt+v l d l ~-d 2 )c 


(17) 


In order to evaluate the running of n it is convenient to 
consider the flow equation for the functional derivatives 
of F fc . In particular we derive the flow equation | 6 ]) with 
respect to 9 and 9 and extract the term proportional to 
the momentum square. The central object in the flow 
equation is the Hessian of the EAA which in field space 
and at zero field reads 


G) 


/ 0 le g 0 0 0 \ 

lee lee 0 0 0 

0 0 7™ 0 0 

0 0 0 0 —Ice 

V 0 0 0 7ec 0 / 


Thus the regularized propagator G k = ^F^ + Rk'j 


~ S 2 T k 1 


G k 


<ST 


( 2 ) 


69 


■G k 


6T 


( 2 ) 


69 


■ G k ■ d f Rk 


+G k ■ <)Tfc 


( 2 ) 


-G fc 


69 

6 2 rf 

6969 


r ^ 


( 2 ) 


69 

G k ■ d t R k 


Gk ■ dtR k 


(18) 


Note that the term in the third line does not contribute 
since there is no such vertex in the ansatz Thus we 
find that 


} a 2 r k 

‘ 6959 


G vv 7^ GggG vv d t Rk,w 


(19) 


where 7 ^ denotes the vertex com ing from the term 
9v l di9 in the action. Equation |l9] ) can be written in 
a diagrammatic fashion as shown in Figure [T] Note that 
there are further contributions to the flow equation for 
the two point function but they vanish upon integration 
over the frequency (see appendix [B] for some details on 
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FIG. 2. Diagrammatic representation of the further contri¬ 
bution to the flow equation of k which vanish after frequency 
integration. 


the frequency integration). These diagrams are shown in 
Figure [2] The final result is: 


Gw')^ Ggg'y^ G vv d t Rk,w — ~ 


dz- 


2 1 

d/ 2-1 


1 27T d / 2 d-l\ Do 

2 


f d+e 


Rk ,1 


(27T) d r(f) d 
2 dtRk,vv ■ 


( 20 ) 


If we use the exponential cutoff ([9]) in the above equation 
we obtain 


Gw'y^ G es "/^ G vv dtRk,m — ~ ^ 


1 f 1 2t r d / 2 d - 1 


( 27 r) d r(|) d 


— k~ e 
2 k 


2^T 

a 


1 + 


d H- s 


( 21 ) 


and for the optimized cutoff (101 


G vv j^G gg ”/ < ' 3 ' , G vv d t Rk,vv — ~ 


1 f 1 2n d / 2 d- 1 

2 V(2-) dr (l) d 

d -\- s 


*k-- 
2 


( 22 ) 


The fact that beta functions coming from (21) and (22) 
are slightly different has to be expected since the beta 
functions are scheme dependent and using different cut¬ 
offs corresponds to a change of scheme in the FRG lan¬ 
guage. Scheme independence has to be recovered in phys¬ 
ical quantities like the critical exponents as we will check 
in a moment. In this case it is convenient to parametrize 
the flow via the coupling g = Dq/k and to introduce its 
dimensionless analogue g = gk~ E . Note that the flow 
equation (20) has the following generic form 

d t (|) = —ADok~ e 


where A is a constant which depend on the dimension 
and on e. In term of the dimensionless coupling g the 
above equation reads: 


d t g = -eg + 2Ag 2 . 


The IR fixed point thus reads g* = e/ (2 A) and the crit¬ 
ical exponent is /3g = e independently of the cutoff. 

Let us compare our procedure with the one adopted in 
[24] . In [24] the authors set up a field theoretic frame¬ 
work analogous to the one we have seen in section [IT] The 
only difference being the fact that the ghost fields are not 
present. The one loop correction to the molecular diffu- 
sivity is computed and the diagram is UV convergent. 
However when e->0a divergence appears. In full anal¬ 
ogy with dimensional regularization the authors perform 
an e —expansion and subtract the poles. This leads to a 
prescription for the RG flow which can be used to obtain 
a resummed result. When one performs this e —expansion 
one notes that the scale m disappears from the leading 
expression proportional to 1/e and disappears from the 
flow equations. In our case we simply set m = 0 from 
the very beginning. 

Note that in comparison to other works in the liter¬ 
ature there are two main differences. First the velocity 
covariance ([ 2 ]) adopted in this work and in [24] differs 
from the one used in other works since the coefficient D 0 
in ([ 2 ]) is replaced with DqE (see for instance EH HEBE!]). 
Moreover the problem can be represented using the Ito 
discretization prescription in place of the Stratonovich 
one. With this definition the velocity covariance is finite 
in the limit e —> 0 and the previously mentioned diver¬ 
gence for e —> 0 is not present [28]. In this case the RG 
flow of the action is trivial and non-trivial flows can be 
found for composite fields like {d6 ) n . We refer to [58] for 
more details on the use of Ito prescription and the RG 
flows found in this case. 


Finally in [21] it is argued that the beta function of 
g is one-loop exact. To understand how this should be 
interpreted in the FRG scheme let us consider again the 
ansatz © and plug this in the flow equation for the two 
point function. We want to understand if further terms 
are generated by the RG flow or not. It is immediate to 
note that no further term is generated: the “trace factor” 
found by contracting the vertices with the propagator is 
— P 2 (l — X 2 ) and the propagator G e § (which depends 
on pi + qi) loses its P —dependence after the frequency 
integration. As a consequence the only P —dependence 
comes from the trace factor. This means that the l.h.s. 
and r.h.s. of (19) are proportional to P 2 , a very uncom¬ 
mon fact in FRG computations. One can be even more 
general and start from the following ansatz 


T k = 


JH* 


1 ^ 


1 


dt + v l di + -/ (-d 2 ) 6 - -0D S 6 + -vD~ l v 


+c (d t + v i d i + -f(-d 2 )j c 


Repeating the same steps the flow equation itself tells 
that / (— d 2 ) = —d 2 . We interpret this as an argument 
telling us that the beta function of g is one-loop ex¬ 
act. Clearly one may consider even more general ansatzs 
which are consistent with all the symmetries we analyzed 
in appendix [A] For example one can add higher powers 
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in the fields via terms like: 

0 2 f (v t , -9 2 ) 0 2 + 20c/ (V t , -9 2 ) 0c. 

We will not consider such terms but these terms should 
be taken into account in a more complete computation. 


V. RG FLOW OF COMPOSITE OPERATORS 

In this section we study the scaling of the structure 
function at the IR fixed point found in section[lV] We will 
first study the renormalization of the operator (9090)" 
and then we will move to consider 0 m . The results found 
in this section will be used in section IVII to obtain the 
correction to the scaling of S 2n . 

A. The operator (d9dd) n 

In section |III| we argued that the scaling of the struc¬ 
ture functions S 2n is encoded in the scaling of the com¬ 
posite operator (9090)". The argument is based on the 
OPE of S‘ 2 n ■ Let us stress that if we take into ac¬ 
count only the operator (9090) we are making an ap¬ 
proximation since other operators also contribute, e.g.: 
(9 2 09 2 0) . In this sense we apply the logic of the deriva¬ 
tive expansion to the composite operator coming from the 
OPE of S 2n - 

From now on we focus on the RG flow of the opera¬ 
tor O = (9090)". To compute its running we introduce 
a coupling Zo which multiplies O and use equation (J7|. 
To extract the running of Zo we perform 2 n functional 
derivatives with respect to 0 and set all the fields to zero. 
Thus l.h.s. is proportional to P 2n and has to be com¬ 
pared with the terms also proportional to P 2n obtained 
by taking 2 n functional derivatives of the r.h.s. of ([7|. 
Note that the Hessian of the operator O appears in equa¬ 
tion (7]) and 2n — 2 functional derivatives will act on it 
while the remaining two will act on the regularized prop¬ 
agators of the EAA. Moreover the l.h.s. of the equation 
is proportional to dtZo while the r.h.s. to Zo so it will be 
easy to extract the anomalous dimension dtZo/Zo■ The 
flow equation obtained in this way can be represented di- 
agrammatically as in Figures [3j [4] and [5] where the black 
dot denotes the insertion of the Hessian of the compos¬ 
ite operator with 2n — 2 functional derivatives attached. 

We limit ourselves to consider the leading term in the 
e —expansion of the structure function which amounts to 
a one loop computation. Our main objective is to test the 
tools introduced so far and reproduce the known leading 
correction. In each diagram the “trace factor” between 
the vertices and the velocity projector gives P 2 (l — A 2 ), 
where A is the cosine between the external and the loop 
momentum. The Hessian of the composite operator reads 
0(2) = 2nP 2(n ~ 1) Q 2 [1 + A 2 2 (ra - 1)], where P and Q 
are the modulus of the external and loop momentum re¬ 
spectively. In each diagram the product of the “trace 



FIG. 3. Diagrammatic representation of a contribution to 
flow equation for composite operators. 



FIG. 4. Diagrammatic representation of a contribution to 
flow equation for composite operators. 


factor” and of 0 (2 ) gives a term of order P 2n . Thus ne¬ 
glecting all the P—dependence in the regularized prop¬ 
agators we obtain the running of Zo and the angular 
integral factors out. The contribution to dtZo/Zo from 



FIG. 5. Diagrammatic representation of a contribution to 
flow equation for composite operators. 








B. The operator 6 m 


the diagram of Figure [3] is 


1 

2 



1 / 27T d / 2 \d~ 1 


<^ d / 2 - 


p 


k,ee 


2n 


d + 2n 
d + 2 

1 


d+e 

2 


+ -Rfc 




Do 

K 

~2 9tRk,vv • 


Furthermore both the diagrams of Figures [4] and [5] give 
the following contribution to dtZojZo-. 


1 

2 


1 ( 2t r d / 2 \ rf-1 

2 VrTiyy 



<tez d / 2 


i±£. . r, 

z 2 + n k ,v 


2 n 


d + 2n 
d + 2 

1 


( z + R k ,ee) 


l£o 

2 K 




k,ee ■ 


The integrals can be solved analytically employing the 
optimized cutoff (101 and, summing all the contributions, 
the leading contribution for small e is 


d t Z 0 1 

f 1 

( 27T d / 2 \ d~ 1 

U n d+2n ] 

Zo 2 

2 

d . 

d + 2 _ 


At the leading order in e we have: 

1 / 2n d D \ 

2Vr(l)J^ 

Therefore we obtain 
dtZo 

Zo 9* 

Although this computation is shown using the optimized 
cutoff we checked that the fixed point anomalous dimen¬ 
sion above is very stable under changes of the cutoff ker¬ 
nel R k . 

Let us note that differently from what we have seen 
in section |IV| the r.lr.s. of the flow equation for com¬ 
posite operators generates new terms. Indeed it is pos¬ 
sible to verify that terms of higher power in P are gen¬ 
erated expanding those propagators which are function 
of (pi + qi). This means that, within our approximation, 
the system is not closed. This is a common problem in the 
FRG approach and the crucial point is to parametrize the 
quantities of interest with sufficiently many terms. Here, 
retaining just one operator, we did a rather crude approx¬ 
imation and this precludes for instance the possibility of 
extracting correct values of the anomalous dimension for 
large values of e. However we stress that hopefully such 
difficulties can be solved by employing larger truncations 
as usually happens in the FRG applications. In order to 
understand if this methodology is efficient for the Kraich- 
nan model one should increase the truncation and check 
the stability of the results. This however is beyond the 
scope of the present work. 


d + 2n 
1 d + 2 



Besides the anomalous dimension of the operator 
( d0d6) n it is useful to consider also the anomalous di¬ 
mension of operators 9 m . It turns out that the anoma¬ 
lous of 9 m dimension is zero. To see this we consider 
the flow equation for the composite operator 9 m which is 
fully analogous to the one for ( 89d9 ) n . One only needs 
to change the “composite operator vertex” in the dia¬ 
grams. This vertex is simply given by m (to — 1) (from 
the Hessian of 9 m ). Now we recall that the trace over 
the vertices and the projectors in these diagrams gives 
an overall factor of P 2 . Since the definition of the run¬ 
ning of Zg™ is taken at zero momentum we see that the 
r.h.s. of the flow equation for Zg™ is zero. 

As an aside comment let us note that a non-trivial 
result is found considering the Kraichnan model for a 
compressible fluid as done in [33] . In this case the action 
reads 


Sg = 


9 (d t 9 + d x {v l 9) - 1 0 2 9 ) - \jD s 9 + \vD~\ 


with 


D„ = 


D n 


(27 


jS (t — t') / d d k- 


aP; 


( k 2 + to 2 ) 2 


l 3 ikx\—ikx2 
+§ 


where Pj :) and PL are the transverse and longitudinal 
projectors respectively. This model has been studied 
in J44j and we have checked at one loop that the fixed 
point and the anomalous dimensions d t Zgm /Zgm are re¬ 
produced by our framework. 


VI. SCALING OF THE STRUCTURE 
FUNCTIONS 


We use the results of sections |IV] and [V] in order to 
derive the running of the structure functions. To achieve 
this we will consider the OPE of [23] . Let us note 
that 


S 2n = ((0 (x)-9(y)) 

2 n 

= E 


2 n\ 


k =0 


2 n 
k 


(9(x) 2n ~ k 9(y) k ). 


Repeating the same steps which lead to equation (16) we 
obtain that all the terms in the above sum satisfy the 
following equation 

[rd r + (2 - e) Td r - n (2 - e) - Md M ] (9 2n ~ k 9 k ) = 0 
where r = \x — y\. This tells us that 

~ r"< 2 - £ ) / (Mr) . 


In order to understand the M —dependence in S 2n we 
consider the OPE. In section III we argued that the lead¬ 
ing term in the OPE of S 2n is given by (d9d9) n thus 

S 2n ~ (C (r) (d989) n ) 
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where C (r) is the first Wilson coefficient of the expan¬ 
sion, a function of (x — y) /2, and ( d6d0) n is a composite 
operator in (x + y) /2. Thus we have 

S 2n ~C(r)((dedO) n ) 


where (( dOdO) n ) is a constant in the sense that it does not 
dependent on any coordinates. Now we consider equation 
(16) for the above correlation and deduce that 


(- MBm + ne + Jq) (( dOdO) n ) = 0 


where we used the fact that (( dddd) n ) does not depend 
on any coordinates. This tells us that 

((d9d9) n ) ~ M ne+1 ° . (23) 


Finally we suppose that the Wilson coefficient C (r) does 
not depend on M in the limit M —>■ 0, i.e.: that C (r) 
is regular in M. Note that the limit M —► 0, i.e.: 
Lp 3> r, is justified by the fact that we are studying 
the inertial regime of the system. This tells us that the 
whole M —dependence in S 2n is given by expression (23). 
Therefore we conclude that 


S- 


2 n 




(Mr) 


ie+7 o 


(24) 


where 


ne + 'y * 0 


2 n (n — 1) 
d + 2 £ 


This is the leading correction found in El® via the zero 
mode approach and in |24l via the e—expansion. 

Let us note that we can derive the scaling of S 2n 
also via the following reasoning. Since (( d6dd) n ) does 
not depend on r we can limit ourselves to study the 
r —dependence in S 2n studying the Wilson coefficient 
C (r). We suppose that C (r) is regular in M and con¬ 
sider the M —> 0 limit. The Callan-Symanzik equation 
for the Wilson coefficient is: 


(/id^ + Pgd g - 70 ) C = 0 

where we used the fact that the anomalous dimension of 
the composite operators 0 2n ~ k and 9 k are zero. Once 
again we repeat the steps which lead to equation (16) 
and obtain 


(■ rd r - 2n + (2 - e) Td r - 7o - Md M ) C = 0. 

Thus we obtain 

C (r) ~ r 2n+1 ° 

where the exponent is exactly the one coming from the 
sum of the exponents in (24). 


VII. SUMMARY AND OUTLOOKS 

We have explored the application of the FRG frame¬ 
work to the Kraichnan model. In particular we studied 
the scaling of the structure functions in the Kraichnan 


model and reproduced the leading corrections by means 
of FRG techniques. From the field theoretical point of 
view the structure functions are correlations of composite 
operators. Thus in section |III| we recalled how to com¬ 
pute the running of composite operators using the flow 
equation for the EAA. In section|TV|andfV|we studied the 
IR fixed point of the model and the anomalous dimen¬ 
sions of various composite operators at the fixed point. 
The scaling of the structure functions in the Kraichnan 
model has been found in section |VT| using the OPE. 

The path integral formulation adopted in this work has 
been explained in section [n] while the associated sym¬ 
metries and some non -renormalization theorems are dis¬ 
cussed in appendix [A] In the path integral formulation 
of the problem a determinant naturally arises and we 
kept track of this by expressing the determinant via the 
integration over some ghost fields. The ghost sector car¬ 
ries two new symmetries which we called BRS symmetry 
(mixing 9 and c) and ghost symmetry. These symmetries 
are very similar to the ones present in the path integral 
formulation for classical mechanics ED] ■ In the latter case 
more symmetries are present, this is due to the fact that 
the equations of motion for the Kraichnan model are not 
Hamiltonian. Let us note that the Navier-Stokes equa¬ 
tion can be formulated as Hamiltonian equations of mo¬ 
tion setting the viscosity to zero and using the so called 
Clebsh variables. It may turn out interesting to study the 
role of the symmetries and of the RG flow in this context. 
Moreover, as far as classical mechanics is concerned, the 
ghost fields carry information on the Lyapunov exponents 
of the system and thus on the chaoticity/stochasticity of 
the system m and we feel that they could play a role 
also in the path integral approach to turbulence. Here 
we limit ourselves to note that in the Kraichnan model 
the equation of motion for the ghost c corresponds to un¬ 
forced solutions which are the ones involved in the zero 
mode approach. We hope to study these topics in the 
future. 

Finally we would like to stress that the techniques 
adopted in this work can be applied to more complicated 
systems, for instance the randomly stirred Navier-Stokes 
equation. On general grounds this will require to fully 
exploit the functional character of the flow equation in 
order to be able to account for an infinite dimensional 
ansatz and not to rely on the derivative expansion. This 
has been achieved successfully in the case of the Kardar- 
Parisi-Zhang equation in m • In the Kraichnan model 
the situation is somehow simpler since the IR fixed point 
is perturbative for small e and the beta function of g 
is one loop exact |2¥|. Nevertheless we believe that the 
logic adopted here, especially with regards to composite 
operators, can be used in other contexts and more work 
is needed in this direction. The advantage of the FRG 
framework is that it does not rely on a perturbative ex¬ 
pansion and that genuinely non-perturbative results can 
be obtained provided that the EAA is parametrized with 
sufficient accuracy. In particular we hope that a simi¬ 
lar framework to the one adopted in this work may be 
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applied to more complicated situations, like the Navier 
Stokes equation. 
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Appendix A: symmetries 

In this appendix we consider the symmetries associated 
to the Krainchnan model. We first consider the sym¬ 
metries associated to the action Sg which can be easily 
extended to the ghost action. Finally we consider sym¬ 
metries coming from non trivial cancellation in Sg + S g h- 
From these symmetries we will be able to derive some 
non-renormalization theorems. Such program has been 
carried out in great detail for the Kardar-Parisi-Zhang 
equation in [21| and more recently for the Navier-Stokes 
equation in m■ In these cases the analysis showed that 
symmetries entail important constraints on the RG flow. 
Here we perform a similar analysis also including the 
ghost sector. 

The action is symmetric under a sort of Z 2 transforma¬ 
tion: 9 —9 and 9 —> —0. This tells us that there is an 

even number of 9 and 9 fields in any monomial of the ef¬ 
fective action. Moreover it is easy to verify the invariance 
under constant shift in 9 

9 —¥ 9 + u. 

Clearly this property is enjoyed also by the effective ac¬ 
tion. In analogy with what has been done in [21] we 
consider a time dependent shift: 

9 —> 9 + u (t) . 

This is not a symmetry but the bare action is non¬ 
invariant just by a term which is linear in field space. 
In particular the only time derivative appears in 

8d t 9 -a 8d t 8 + 8d t u ( t ) . 

Since the non-invariance of the action is given just by a 
term which is linear in the held we can derive the follow¬ 
ing Ward Idendity (WI)Q 

0 = (J ■ u (t) — J 9d t u (t)) 

= ■«(*)-/ 8d t u{t)j). 

1 We use the following conventions: exp W = f exp [—S + J ■ <fi] 
and T = f Jip- W. Thus J = S V F. 


This equation tells us that the term in the bare action 
f 9dt9 does not get renormalized. Moreover the action 
and the effective action also enjoy Galilei invariance as 
we will see in a moment. Thus we can infer that the term 
/ 9 (pt + v l di ) 9 is not renormalized. 

We also consider a time dependent shift in the response 
held 9 —> 9 + e(t) and we have 

0 = (- J e(t) (d t 9) + e ( t ) (ff)9 + J g e t ) 

~%' £t = ~J ( at6l ) +£ W 

This entails the non-renormalization of the term 

Now we consider the Galilei invariance of the model, 
the transformation is: 

v (t, x) —> v' (t, x') + c = v' ( t , x — ct) + c 
x —► x' = x — ct 
t -A t' 

9 (t,x) —> 9 (t,x') 

9 (t, x) —> 9 (■ t , x') . 

This invariance tells us that time derivative actually ap¬ 
pears via Galilei covariant derivative of the form dt+v l di. 
Note that the term quadratic in the velocities is not 
straightforwardly invariant. This is due to the fact that 
the average velocity is set to zero, in order to see directly 
Galilei invariance we must express the kinetic term as 
(v — (v)) D v (v — (v)). More precisely it turns out very 
useful also to consider the following time gauged Galilei 
transformation: 

v (■ t , x) —> v' (t , x') + c(t) = v' (t,x — c(t)) + c ( t) 
x —> x' = x — c (t) 
t —y t 

9 (t,x) —> 9{t,x') 

9 (■ t , x) 9 (■ t , x'). 

We observe that 

(dt + vdx) 9 —> [dt9 t x — c{t ) “1“ d x 0t,X—c(t) ‘ ( d t c)\ 

“b (Wz-c(t) T (^)) ^x9t,x—c{t ) 

= [(d t + vd x ) 9}'. 

Note that this is a symmetry of the equation ([I]) itself 
which is not preserved by the averaging over the veloc¬ 
ities. Indeed performing a transformation in the path 
integral we obtain 

0 = (— J vD v Sv + J v Sv). 

This entails the non-renormalization of the “kinetic” 
term for the velocity fields. Thus we observe that the ef¬ 
fective action is invariant under the time-gauged version 
of the Galilei transformations except for the quadratic 
term in the velocity field. All these non-renormalization 
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theorems are known in the perturbative framework, see 

M- 

So far we have been considering symmetries which do 
not rely on the presence of the ghost. Indeed Sg and S g h 
are both simultaneously invariant. Now we examine sym¬ 
metries which exploit a non-trivial cancellation between 
Sg and Sgh- Let us note that the action can be rewritten 
as 

S = - ^ ODgO + OO v 6 + cO v c + l>vD~ 1 11 v 

where O v = d t + v l di — |<9 2 . It is easy to check that 
0 — > 6 + ec, c — ^ c — e9 is a symmetry which we call 
BRS symmetry in analogy to nomenclature in [3D]. This 
symmetry has been discussed in [36j and we briefly re¬ 
view it here. The symmetry transformation is linear and 
thus the BRS invariance is present also in the effective 
action. To derive some useful results we can consider 
(SbrsO) = 0 for some particular operators. For in¬ 
stance we can consider (5 ( cO)) = {(—00 + cc)) =0 and 
(6 (cd)) ~ (66) = 0. Finally we have also the following 
symmetry: c —> c — ec and c —> c + ec. This symmetry 
entails the ghost number conservation. 


Appendix B: some frequency integrals 

In this appendix we provide a list of the typical fre¬ 
quency integral we have encountered. These integral 
come from the multiplication of various 60 —propagator 
entering in each diagram. The regularized propagator 


has the following form: 

1 _ + V 

±iLO + V ~ u 2 + V 2 ‘ 

First of all we consider 

G(q) = J du (—iu + vq 2 ) 1 . 

This integral is divergent but its principal value is well 
defined: 



J dco (—iu> + V) 1 




iuj + V 
u 2 + V 2 


Moreover we need 


1 

27r 


1 

27T 


( —-— 


2 


= 0 


dui 


1 1 1 
ioo + V —ioo + V iuj + V 


1 1 

4V 2 ' 


The first integral above is somewhat unusual in the this 
type of computations, since one typically finds 

1 f 1 1 _ 1 

27t J —iuj + V iuj + V 2V 
This happens because our diagrams contains not only 
60 —propagator but also velocity propagators (which are 
delta in time). 
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